200 A GeV Au+Au collisions serve a nearly perfect quark-gluon liquid 
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A new robust method to extract the specific shear viscosity (r?/s) qgp of a Quark-Gluon-Plasma 
(QGP) at temperatures T c < T < 2T C from the centrality dependence of the eccentricity-scaled elliptic 
flow v^je measured in ultra-relativistic heavy- ion collisions is presented. Coupling viscous fluid 
dynamics for the QGP with a microscopic transport model for hadronic freeze-out we find for 
200 A GeV Au+Au collisions that v^jz is a universal function of multiplicity density (l/S)(dN c h/dy) 
that depends only on the viscosity but not on the model used for computing the initial fireball 
eccentricity e. Comparing with measurements we find 1 < 47t(7?/s)qgp < 2.5 where the uncertainty 
range is dominated by model uncertainties for the values of e used to normalize the measured «2- 
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Ever since heavy-ion collision experiments at the Rcl- 
ativistic Heavy-Ion Collider (RHIC) demonstrated the 
creation of color-deconfined Quark-Gluon Plasma (QGP) 
exhibiting almost ideal fluid dynamical collective behav- 
ior [1-3], with viscosity per entropy density rj/s ap- 
proaching the KSS lower bound 7^^: [4, 5], an ac- 
curate extraction of the QGP transport coefficients, es- 
pecially its shear viscosity (?7/s)qgp, from experimental 
measurements has been of great interest [6-8]. A small 
•q/s is generally considered to be evidence for the on- 
set of a strongly-coupled deconfined plasma early in the 
evolution of the collision. RHIC is the first accelerator 
to provide sufficient beam energy for the QGP to live 
long enough for flow observablcs to become sensitive to 
its intrinsic transport properties. Simulations based on 
both viscous fluid dynamics and quark-gluon transport 
theory [6-8] have established that the elliptic flow gen- 
erated in non-central heavy-ion collisions is particularly 
sensitive to the shear viscosity r\j s of the medium. How- 
ever, a quantitative extraction of (?7/s)qgp from elliptic 
flow data requires not only accurate elliptic flow mea- 
surements but also a precise knowledge of the theoreti- 
cal baseline corresponding to zero QGP viscosity. The 
latter, in turn, requires good control over the fluid's col- 
lective response to anisotropic pressure gradients, and a 
realistic microscopic description of chemical and kinetic 
freeze-out during the hadronic stage [9] . A purely hydro- 
dynamic approach that treats both the dense early QGP 
and dilute late hadron resonance gas phases as viscous 
fluids not only requires the introduction of two additional 
parameters, the chemical and kinetic freeze-out tempera- 
tures, which must be separately adjusted to experimental 
data, but ultimately fails [10] because viscous corrections 
due to hadronic dissipation are large [11] and invalidate 
a fluid dynamical approach even if it properly accounts 
for chemical decoupling before kinetic freeze-out [12, 13] 
and for a strong growth [14] of the specific shear viscosity 



77/s in the hadronic stage [15]. 

We here use a newly developed hybrid code (see [10] for 
details) that couples the relativistic (2+l)-dimensional 
viscous fluid algorithm VISH2+1 [8] to the microscopic 
hadronic scattering cascade UrQMD [16] via a Monte Carlo 
interface [17]. For the QGP fluid we assume constant 
77/s for T C <T<2T C [18]. We switch from a hydrody- 
namic description of the QGP to UrQMD at tempera- 
ture T sw = 165 MeV, adjusted to reproduce the chemi- 
cal freeze-out temperature measured in RHIC collisions 
[19] and the highest T for which we have a valid micro- 
scopic description. By giving us full microscopic control, 
without additional parameters, over the complex hadron 
kinetic freeze-out our hybrid model opens the door for 
quantitatively exploring the transport properties of the 
earlier QGP phase using measured final hadron spectra. 

For the hydrodynamic evolution above T sw we use the 
state-of-the-art equation of state (EOS) s95p-PCE based 
on recent lattice QCD results [13]. The remaining model 
uncertainties arise mainly from the initial conditions of 
the hydrodynamic evolution, including the starting time 
To and initial transverse flow velocity. While these can- 
not be directly measured and require model input, they 
are tightly constrained by experimental information on 
the final state [1, 20]. Modeling the QGP as an ideal 
fluid with r]/s = and zero initial transverse flow requires 
an early start at To =0.4fm/c. Non-zero shear viscosity 
adds to the transverse pressure [6-8] , generating stronger 
radial flow. The same final flow can then be reached 
with later starting times, giving the system more time 
for thermalization. We find [21] that the shapes of the 
measured pion and proton pr-spectra are well reproduced 
with the following parameter pairs (77/s, ctq): (0, 0.4 fm), 
(0.08,0.6fm), (0.16,0.9fm), and (0.24, 1.2 fm). 

For each choice of To, the initial energy density is renor- 
malized to yield the same final charged hadron multiplic- 
ity dNch/dy in central Au+Au collisions. Its distribution 
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FIG. 1. (Color online) (a) Centrality dependence of the charged hadron rapidity density per participant pair 
(dN ch /dy)/(N ps , It /2). Experimental data are from STAR [33] and PHOBOS [37], using dN ch /dy = 1.16 dN ch /dri for PHO- 
BOS. Theoretical lines are explained in the text, (b) Eccentricity-scaled elliptic flow vi/e as function of multiplicity density 
(1/S) dN c h/dy, for different values of (t)/s)qgp- Here and in Fig. 2 V2 is integrated with the same cuts as in the STAR data 
[38]: 0.15GeV/c<pr <2GeV/c, \q\ < 1. The overlap area S is always from the same initial state model as the eccentricity e 
(see text). Note the universality of this theoretical relation, independent of the model used for calculating e and S. Panels (a) 
and (b) use the same colors and symbols but for clarity not all corresponding curves are shown in both panels. 



in the transverse plane is determined (via the EOS) from 
the initial entropy density distribution s(r,To;b) which 
we compute, alternatively, from two geometric models 
discussed below. For the shear viscous pressure tensor 
we use Navier-Stokes initial conditions [8], noting that 
the system loses memory after a few relaxation times t w 
where = f£ = 0(0.2 fm/c) [22]. We ignore bulk vis- 
cosity due to its small effect on p^-spectra and vi [23] . 

The key driver for the elliptic flow generated in the col- 
lision is the initial source eccentricity e = where x 
and y label the coordinates along the short and long ma- 
jor axes of the fireball in the transverse plane, e is com- 
puted from the initial entropy density after thermaliza- 
tion [24] . For a quantitative comparison with experiment 
we account for event-by-event fluctuations of e [25] as 
follows: For each impact parameter, we generate an en- 
semble of initial entropy density distributions by Monte 
Carlo (MC) sampling an analytic model of the collision 
geometry, recentering and rotating each distribution such 
that its short major axis x aligns with the direction of 
the impact parameter. The plane defined by the short 
major axis and the beam direction (x— z plane) is called 
"participant plane" , and the eccentricity using this def- 
inition of x is denoted as e pa rt- Superimposing many 
such events yields a relatively smooth input distribution 
for hydrodynamic evolution, with an average eccentricity 
(e part ). The resulting elliptic flow is interpreted as the 
event-average (v 2 ) for the selected centrality class. 

Experimental methods for extracting the elliptic flow 
[26] typically do not yield (1*2) ■ For example, the 2- 
particle cumulant, denoted by «2{2}, includes event-by- 
event flow fluctuations, plus so-called "non-flow" contri- 



butions that are outside the purview of hydrodynamics 
[27, 28]. Fortunately, recent work [28] removed these fluc- 
tuation and non-flow contributions from the measured el- 
liptic flow, thereby providing experimental values for (i^) 
that can be normalized by (£ pa rt) for a direct comparison 
with theory. In the absence of non-flow, t>2{2} ~ y (uf); 

assuming [28] VT^IT - ?i^Ty\/ ( £ pa.rt> [29], the experi- 
mentally determined left side can then again be compared 
with the theoretically computed right side. We show such 
a comparison below to check consistency. 

To compute the initial entropy density distribution in 
the transverse plane we use MC versions of the Glauber 
[31] and fKLN [32] models; for a detailed description 
of our procedure sec [24]. The models are tuned to 
reproduce the measured collision centrality dependence 
of dN^/dy. Figure 1(a) shows that, for all permissible 
combinations of To and 77/s and both MC-Glauber and 
MC-KLN models for the initial density distribution, the 
measured centrality dependence of dN c h/dy is well repro- 
duced. The same holds for the slopes of pion and proton 
spectra at all centralities [21]. (Following STAR [33], 
dN^/dy does not include charged hyperons and weak 
decay products.) Two additional curves for initial MC- 
Glauber and MC-KLN densities with uniformly reduced 
(by ~ 10%) final multiplicities arc shown to demonstrate 
that, as long as the overall trend is preserved, small dif- 
ferences in dNch/dy extracted from STAR, PHOBOS and 
PHENIX measurements do not influence our conclusions. 

Figure 1(b) shows the key theoretical result of the 
present study: the relation between eccentricity scaled el- 
liptic flow v-i/e and multiplicity density (l/S)dN c h/dy is 
approximately universal (at least for fixed y/s), depend- 
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FIG. 2. (Color online) Comparison of the universal V2(rj/s)/e vs. (l/S)(dN c h/dy) curves from Fig. 1(b) with experimental 
data for {V2) [28], ^{2} [38], and dN^/dy [33] from the STAR Collaboration. The experimental data used in (a) and (b) are 
identical, but the normalisation factors (e pa rt) and S used on the vertical and horizontal axes, as well as the factor (£p art ) 1//2 
used to normalize the «2{2} data, are taken from the MC-KLN model in (a) and from the MC-Glauber model in (b). Theoretical 
curves are from simulations with MC-KLN initial conditions in (a) and with MC-Glauber initial conditions in (b). 



ing only on the value of 77/s for the QGP but not on any 
details of the model from which s and S = nyJ (x 2 )(y 2 ) 
are computed. To good approximation, switching be- 
tween initial state models shifts points for a given colli- 
sion centrality along these universal curves, but not off 
the lines. For example, reducing the final multiplicity 
by renormalizing the initial entropy density shifts the 
points towards the left but also downward because less 
elliptic flow is created, due to earlier hadronization. The 
significantly larger (e pa rt) from the KLN model gener- 
ates more v 2 than for the Glauber model, but the ratio 
v-ije is almost unchanged. Slightly larger overlap areas 5* 
for the KLN sources decrease (l/S)(dN c h/dy), but this 
also decreases the initial entropy density and thus the 
QGP lifetime, reducing the ratio V2/S; the result is a 
simultaneous shift left and downward. Early flow [34] 
(r =0.4fm/c for i]/s = 0.08) increases v 2 /s by ~ 5%, 
but the separation between curves corresponding to ry/s 
differing by integer multiples of l/(47r) is much larger. 
Only in very peripheral collisions is the universality of 
v 2 /e vs. (l/S)(dN ch /dy) slightly broken [35]. 

The clear separation and approximate model- 
independence of the curves in Fig. 1(b) corresponding 
to different (?7/s)qgp values suggests that one should be 
able to extract this parameter from experimental data. 
However, only v 2 and dN^/dy are experimentally mea- 
sured whereas the normalization factors e and S must be 
taken from a model. Figure 2 shows a comparison of the 
theoretical curves from Fig. 1(b) with STAR data nor- 
malized by eccentricities and overlap areas taken from 
different initial state models that were all tuned to cor- 
rectly reproduce the centrality dependence of dN^/dy 
shown in Fig. 1(a) [39]. Since, for the same model, 
the eccentricities and overlap areas depend somewhat on 



whether they are calculated from the initial energy or en- 
tropy density, the same definitions must be used in theory 
and when normalizing the experimental data. 

Both panels of Fig. 2 show the same data, in panel 

(a) normalized by e, S from the MC-KLN model and in 

(b) with the corresponding values from the MC-Glauber 
model. The theoretical curves arc from the same models 
as used to normalize the data. The figure shows that 
comparing apples to apples matters: When comparing 
the data for v 2 {2} / {e^ avt ) 1/2 with those for («2)/(£ P art), 
the former are seen to lie above the latter, showing that 
non-flow contributions (which cannot be simulated hy- 
drodynamically) either make a significant contribution 
to v 2 {2} or were overcorrected in (v 2 ) [28], especially in 
peripheral collisions. The extraction of rj/s from a com- 
parison with hydrodynamics thus requires careful treat- 
ment of both fluctuation and non-flow effects. 

The main insight provided by Fig. 2 is that the the- 
oretical curves successfully describe the measured cen- 
trality dependence of v 2 /e, i.e. its slope MS cl function 
of dN^/dy, irrespective of whether the measured elliptic 
flow is generated by an initial MC-KLN or MC-Glauber 
distribution. To the best of our knowledge, the hybrid 
model used here to describe the dynamical evolution of 
the collision fireball is the first model to achieve this. The 
magnitude of the source eccentricity (and, to a lesser ex- 
tent, of the overlap area) disagrees between these two 
models, and this is the main source of uncertainty for 
the value for (??/s)qgp extracted from Fig. 2. Both the 
Glauber and KLN models come in different flavors, de- 
pending on whether the models are used to generate the 
initial entropy or energy density. We have checked that 
the versions studied here produce the largest difference 
in source eccentricity between the models. In this sense 
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we are confident that Figs. 2(a) and (b) span the realistic 
range of model uncertainties for e and S. 

We conclude that the QGP shear viscosity for 
T c <T<2T C lies within the range 1 < Att(?]/s)q GP < 2.5, 
with the remaining uncertainty dominated by insufficient 
theoretical control over the initial source eccentricity e. 
While this range roughly agrees with the one extracted 
in [7] , the width of the uncertainty band has been solid- 
ified by using a more sophisticated dynamical evolution 
model which eliminates most possible sources of error 
that the earlier analysis [7] was unable to address. Small 
bulk viscous effects [23] and proper event- by- event hy- 
drodynamical evolution of fluctuating initial conditions 
[30] may slightly reduce the ideal fluid dynamical base- 
line, while pre-equilibrium flow may slightly increase it. 
Although this should be studied in more quantitative de- 
tail, we expect the quoted unccrtainy band for (t?/s)qgp 
to shift, after cancellations, by only a few percent. 
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